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MAGNETIC WELLS IN DIMENSION THREE 


B. HELFFER, Y. KORDYUKOV, N. RAYMOND, AND S. VU NGOC 


Abstract. This paper deals with semiclassical asymptotics of the three- 
dimensional magnetic Laplacian in presence of magnetic confinement. Us¬ 
ing generic assumptions on the geometry of the confinement, we exhibit 
three semiclassical scales and their corresponding effective quantum Hamil¬ 
tonians, by means of three microlocal normal forms a la Birkhoff. As a 
consequence, when the magnetic field admits a unique and non degenerate 
minimum, we are able to reduce the spectral analysis of the low-lying eigen¬ 
values to a one-dimensional ft-pseudo-differential operator whose Weyl’s 
symbol admits an asymptotic expansion in powers of h? . 


1. Introduction 

1.1. Motivation and context. The analysis of the magnetic Laplacian 
(—iKV — A) 2 in the semiclassical limit h —> 0 has been the object of many 
developments in the last twenty years. The existence of discrete spectrum 
for this operator, together with the analysis of the eigenvalues, is related to 
the notion of “magnetic bottle”, or quantum confinement by a pure magnetic 
field, and has important applications in physics. Moreover, motivated by 
investigations of the third critical field in Ginzburg-Landau theory for super¬ 
conductivity, there has been a great attention focused on estimates of the 
lowest eigenvalue. In the last decade, it appears that the spectral analysis 
of the magnetic Laplacian has acquired a life on its own. For a story and 
discussions about the subject, the reader is referred to the recent reviews 
[11,14,24], 

In contrast to the wealth of studies exploring the semiclassical approxi¬ 
mations of the Schrodinger operator —h 2 A + V, the classical picture associ¬ 
ated with the Hamiltonian \\p — A(g)|| 2 has almost never been investigated 
to describe the semiclassical bound states ( i.e . the eigenfunctions of low 
energy) of the magnetic Laplacian. The paper by Raymond and Vu Ngqc 
[25] is to our knowledge the first rigorous work in this direction. In that 
paper, which deals with the two-dimensional case, the notion of magnetic 
drift, well known to physicists, is cast in a symplectic framework, and using 
a semiclassical Birkhoff normal form (see for instance [27, 5, 28]) it becomes 
possible to describe all the eigenvalues of order 0(h). Independently, the as¬ 
ymptotic expansion of a smaller set of eigenvalues was established in [12, 15] 
through different methods which act directly on the quantum side: explicit 
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unitary transforms and a Grushin like reduction are used to reduce the two- 
dimensional operator to an effective one-dinrensional operator. 

The three-dimensional case happens to be much harder. The only known 
results in this case that provide a full asymptotic expansion of a given eigen¬ 
value concern toy models where the confinement is obtained by a boundary 
carrying a Neumann condition on an half space in [23] or on a wedge in [22]. 
In the case of smooth confinement without boundary, a construction of quasi¬ 
modes by Helffer and Kordyukov in [13] suggests what the expansions of the 
low lying eigenvalues could be. But, as was expected by Colin de Verdiere 
in his list of open questions in [7], extending the synrplectic and nricrolocal 
techniques to the three-dimensional case contains an intrinsic difficulty in the 
fact that the synrplectic form cannot be nondegenerate on the characteristic 
hypersurface. The goal of our paper is to answer this question by fully car¬ 
rying out this strategy. After averaging the cyclotron motion, the effect of 
the degeneracy of the synrplectic form can be observed on the fact that the 
reduced operator is only partially elliptic. Hence, the key ingredient will be 
a separation of scales via the introduction of a new semiclassical parameter 
for only one part of the variables. These semiclassical scales are reminiscent 
of the three scales that have been exhibited in the classical picture in the 
large field limit, see [2, 6 ]. They are also related to the Born-Oppenheinrer 
type of approximation in quantum mechanics (see for instance [4, 20]). In 
fact, in a partially semiclassical context and under generic assumptions, a 
full asymptotic expansion of the first magnetic eigenvalues (and the corre¬ 
sponding WKB expansions) has been recently established in any dimension 
in the paper by Bonnaillie-Noel-Herau-Raynrond [3]. 

1.2. Magnetic geometry. Let us now describe the geometry of the prob¬ 
lem. The configuration space is 

R 3 = {q 1 ei+q 2 e 2 +q3e 3 , qj E M, j = 1,2,3}, 

where (e ^)^®! 2 3 is the canonical basis of M 3 . The phase space is 

M 6 = {(q,p) E M 3 x M 3 } 

and we endow it with the canonical 2 -fornr 

(1.1) w 0 = dpi A dqi + dp 2 A dq 2 + dp 3 A dq 3 . 

We will use the standard Euclidean scalar product (•, •) on M 3 and || • || the 
associated norm. In particular, we can rewrite ujq as 

uj 0 ({ui,u 2 ), (vi,v 2 )) = {vi,u 2 ) - (v 2 ,ui), \/ui,u 2 ,vi,v 2 E M 3 . 

The main object of this paper is the magnetic Hamiltonian, defined for all 
(q,p) € M 6 by 

(1-2) H(q,p) = \\p-A(q)\\\ 

where A E C°°(M 3 ,M 3 ). 
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Let us now introduce the magnetic field. The vector field A = (Ai, A 2 , A 3 ) 
is associated (via the Euclidean structure) with the following 1-form 

a = Aidqi + A 2 d <? 2 + A 3 dq 3 

and its exterior derivative is a 2 -form, called magnetic 2 -form and expressed 
as 

da = ( diA 2 -d 2 Ai)dqiAdq 2 +(diA 3 -d 3 Ai)dqiAdq 3 +(d 2 A 3 -d 3 A 2 )dq 2 Adq 3 . 
The form da may be identified with a vector field. If we let: 

B = V x A = (d 2 A 3 - <9 3 A 2 , d 3 A 1 - 5iA 3 , 5iA 2 - d 2 A 1 ) = (B u B 2 , B 3 ), 
then, we can write 

(1.3) da = B 3 dqi A dq 2 - B 2 dq\ A dq 3 + B\dq 2 A dq 3 . 

The vector field B is called the magnetic field. Let us notice that we can 
express the 2 -form da thanks to the magnetic matrix 

/ 0 B 3 -B 2 \ 

Mb = j — B 3 0 Bi 

\ B 2 -B\ 0 J 

Indeed we have 

(1.4) da(U, V ) = (U, M b V) = (U,Vx B) = [U, V, B], V(C7, V) G M 3 x M 3 , 

where is the canonical mixed product on M 3 . We note that B belongs 

to the kernels of A/b and da. 

An important role will be played by the characteristic hypersurface 

X = H~ 1 ( 0), 

which is the submanifold defined by the parametrization: 

K 3 3giA j(q) := (q, A (q)) G M 3 X M 3 . 

We may notice the relation between £, the symplectic structure and the 
magnetic field in the following relation 

(1.5) j*u 0 = da, 
where da is defined in (1.3). 

1.3. Confinement assumptions and discrete spectrum. This paper is 
devoted to the semiclassical analysis of the discrete spectrum of the magnetic 
Laplacian Cti,A '■= (— iKV q — A (q)) 2 , which is the semiclassical Weyl quanti¬ 
zation of H (see (2.1)). This means that we will consider that h belongs to 
( 0 , ho) with ho small enough. 

Let us recall the assumptions under which discrete spectrum actually exist. 
In two dimensions, with a non vanishing magnetic field, a standard estimate 
(see [ 1 , 8 ]) gives 

h [ \B(q)\\u(q)\ 2 dq ^ (C h , A u \ u) , V« G C 0 °°(M 2 ) . 

Jr 2 
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Except in special cases when some components of the magnetic field have 
constant sign, this is no more the case in higher dimension (see [10]). We 
should impose a control of the oscillations of B at infinity. Under this con¬ 
dition, we get a similar estimate at the price of a small loss. This kind of 
estimate actually follows from an analysis developed in [16]. Let us define 

b(q) := ||B(g)||. 

Let us now state the confining assumptions under which we will constantly 
work in this paper. 

Assumption 1.1. We consider the case of M 3 and assume 

(1.6) b(q ) ^ b 0 := inf b(q) > 0, 

<jGK 3 

and the existence of a constant C > 0 such that 

(1.7) ||VB( g )|| <C(1 + &(</)), Vg E M 3 . 

Under Assumption 1.1, it is proven in [16, Theorem 3.1] that there exist 
ho > 0 and Co > 0 such that, for all h E (0, ho), 

(1.8) h{l — Coh±) [ b(q)\u(q)\ 2 dq ^ {C-h.AU \ u) , Mu E C“(M 3 ). 

J R 3 

As a corollary, using Persson’s theorem (see [21]), we obtain that the bot¬ 
tom of the essential spectrum is asymptotically above hb\ , where 

b\ := liminf b(q). 

|g|—>+oo 

More precisely, under Assumption 1.1, there exist ho > 0 and Co > 0 such 
that, for all h E (0, ho), 

(1.9) s eS s(£ft,A) c [^&i(l “ C 0 /j2),+oo). 

Assumption 1.2. We assume that 

(1.10) 0 < 6 0 < &i. 

Moreover we will assume that there exists a point qo E M 3 and e > 0, /3o E 
(bo, b\) such that 

(1.11) {b(q) ^ fa} C D(qo,e), 

where D(qo,e ) is the Euclidean ball centered at qo and of radius e. For the 
rest of the article we let f3o E (&o>A))- Without loss of generality, we can 
assume that qo = 0 and that A(0) = 0 (which can be obtained with a change 
of gauge). 

Note that Assumption 1.2 implies that the minimal value of b is attained 
inside D(qo,s). 

All along this paper, we will strengthen the assumptions on the nature 
of the point qo. At some stage of our investigation, qo will be the unique 
minimum of b. Note in particular that (1.11) is satisfied as soon as b admits 
a unique and non degenerate minimum. 
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1.4. Informal description of the results. Let us now informally walk 
through the main results of this paper. We will assume (as precisely formu¬ 
lated in (l.lO)-(l.ll)) that the magnetic field does not vanish and is confin¬ 
ing. 

Of course, for eigenvalues of order 0(K), the corresponding eigenfunctions 
are microlocalized in the semi-classical sense near the characteristic man¬ 
ifold £ (see for instance [26, 31]). Moreover the confinement assumption 
implies that the eigenfunctions of £} L a associated with eigenvalues less that 
/3oh enjoy localization estimates a la Agmon. Therefore we will be reduced to 
investigate the magnetic geometry locally in space near a point qo = 0 G M 3 
belonging to the confinement region and which, for notational simplicity, we 
may assume to be the origin. 

Then, in a neighborhood of (0, A(0)) G £, there exist symplectic coordi¬ 
nates (xi, £i, X 2 , £ 2 , x 3 , £ 3 ) such that £ = {x\ = £1 = £3 = 0} and (0, A(0)) 
has coordinates 0 G M 6 . Hence £ is parametrized by (x 2 , £, 2 , x 3 ). 

1.4.1. First Birkhoff form. In these coordinates suited for the magnetic ge¬ 
ometry, it is possible to perform a semiclassical Birkhoff normal form and 
microlocally unitarily conjugate £h,a to a first normal form Mn = Op™ (N%) 
with an operator valued symbol 1W depending on (X2, £21X3, £3) in the form 

N h = e 3 +b(x 2 ,&,x 3 )i h +nn,i h , x 2 ,6, *3,6)+ omr, i&n. 

where X/, = h?D 2 i + x\ is the first encountered harmonic oscillator and where 
(H,I,x 2 ,&,X 3 ,&) f*(h,I,X 2 ,&,x 3,^3) satisfies, for I G (0, To), 

\f*(h,I,x 2 ,b,x 3 ,&)\ < C (|/|f + |£ 3 | 3 + ^) • 

Since we wish to describe the spectrum in a spectral window containing at 
least the lowest eigenvalues, we are led to replace X^ by its lowest eigenvalue h 
and thus, we are reduced to the two-dimensional pseudo-differential operator 

A = Op™ where 

N£ ] =e3+Hx 2 ,(;2,x 3 )h + r(h,h,x2,&,x 3 ,z 3 ) + o(h°°,\z 3 n. 

1.4.2. Second Birkhoff form. If we want to continue the normalization, we 
shall assume a new non-degeneracy condition (the first one was the positivity 
of h). 

Now we assume that, for any (x 2 , £2) i n a neighborhood of (0, 0), the func¬ 
tion x 3 1 —> b(x 2 , £ 2 , x 3 ) admits a unique and non-degenerate minimum de¬ 
noted by s(x 2, £2)- Then, by using a new symplectic transformation in order 
to center the analysis at the partial minimum s(x 2 , £ 2 ); we get a new operator 
A/1. 1 ^ whose Weyl symbol is in the form 

4 1 = ^ 2 (x2,6)(^3 + h x l) + hb(x 2 ,&, s(x 2,6)) + remainders, 


with 

( 1 . 12 ) 


^(*2, 62) = (|9|6( x 2,6, s ( x 2 ,6))) 1/4 
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and where the remainders have been properly normalized to be at least for¬ 
mal perturbations of the second harmonic oscillator + hx 3. Since the 
frequency of this oscillator is h~ 2 in the classical picture, we are naturally 
led to introduce the new semiclassical parameter 

h = f& 

and the new impulsion 

C = 

so that 

Op n (£3 + ^ x i) = Op^ (Cl + • 

We therefore get the /i-symbol of A/~^: 

N W _ h 2 v 2 (x 2 , hf 2) (is + X3) + h 2 b(x 2 , h£ 2 , s(x 2, /i| 2 )) + remainders. 

We can again perform a Birkhoff analysis in the space of formal series given 
by S' = J^[x3, 6, h] where ^ is a space of symbols in the form c(h, x 2 , h£ 2 ). 
We get the new operator DJth = Op™ (M/J, with 

= h 2 b(x 2 ,hl 2 ,s(x 2 ,h£ 2 )) + h 2 J h Op™ v 2 (x 2 ,hi 2 ) + h 2 g*(h,J h ,x 2 ,h£ 2 ) 

+ remainders, 

where Jh = Op^ and g*(h , J, x 2 ,6) is of order three with respect 

to (J2, /12). Motivated again by the perspective of describing the low lying 
eigenvalues, we replace Jh by h and rewrite the symbol with the old semi¬ 
classical parameter h to get the operator = Op^ = Op™ , 

with 

U-13) 

= hb(x 2 J 2 ,s(x 2 J 2 )) +h^v 2 (x 2,6) + hg*(hjhjx 2,6) + remainders. 

1.4.3. Third Birkhoff form. The last generic assumption is the uniqueness 
and non-degeneracy of the minimum of the new “principal” symbol 

(^2,6) ^ b(x 2 J 2 ,s(x 2 J 2 )) 

that implies that b admits a unique and non-degenerate minimum at (0, 0, 0). 
Up to an /12-dependent translation in the phase space and a rotation, we are 
essentially reduced to a standard Birkhoff normal form with respect to the 
third harmonic oscillator JCh = h 2 D 2 2 + x\. 

Note that all our normal forms may be used to describe the classical dy¬ 
namics of a charged particle in a confining magnetic field (see Figure 1). 
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Figure 1. The dashed line represents the integral curve 
of the confining magnetic field B = curlA through qq = 
(0.5, 0.6, 0.7) for B (x,y,z) = \/l + x 2 ^j and the full 

line represents the projection in the q-space of the Hamil¬ 
tonian trajectory with initial condition ( qo,Po ) (with po = 
(—0.6,0.01,0.2)) ending at (qi,Pi). 


1.4.4. Microlocalization. Of course, at each step, we will have to provide ac¬ 
curate microlocal estimates of the eigenfunctions of the different operators 
to get a good control of the different remainders. In a first approximation, 
we will get localizations at the following scales £i,£i ,£3 h s {5 > 0 is 
small enough) and X 2 ,^ 2^3 ~ 1- In a second approximation, we will get 
h 5 . In the final step, we will refine the localization by X 2 , £2 ~ h s . 

1.5. A semiclassical eigenvalue estimate. Let us already state one of the 
consequences of our investigation. It will follow from the third normal form 
that we have a complete description of the spectrum below the threshold 
boh + 3u 2 (0, 0)/P. This description is reminiscent of the results a la Bohr- 
Sommerfeld of [17] and [18, Appendix B] (see also [15, Remark 1.4]) obtained 
in the case of one dimensional semiclassical operators. 
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Theorem 1.3. Assume that b admits a unique and non degenerate minimum 
at qo. Denote 


(1.14) 


Hessg 0 6 (B, B) Q j detHess go 6 
262 ’ = y Hess go 6 (B, B)' 


There exists a function k * £ C“(M 2 ) with arbitrarily small compact support, 

1 3 

and k*{h? , Z) = 0{{h + \Z\) 2 ) when (h, Z) —> (0,0), such that the following 
holds. 

1 3 

For all c £ (0, 3), the spectrum of Cn.A below boh+cazh? coincides modulo 
0{h°°) with the spectrum of the operator J-% acting on L 2 (M X ) given by 


Tn = 6 0 ft + era - —h 2 


+ h 


e -K h + k*{H*,K h ) 


with some constant £. 


IC h — h 2 D 2 + x 2 , 


Remark 1.4. The constant f in Theorem 1.3 is given by the formula 

C = l|V^ 2 (0,0)|| 2 , 

where the function u is given in (1.12). Observe also that a = z^ 4 (0, 0). 


Corollary 1.5. Under the hypothesis of Theorem 1.3, let (A m{h)) m>1 be the 
non decreasing sequence of the eigenvalues of Cf, a- For any c £ (0,3), let 

N/, c := {m £ N*; X rn (h) ^ hbo + ca^h 1 *}. 

Then the cardinal of Nj-i,c °f order h~ 2 , and there exist v\, i >2 £ M and 
ho > 0 such that 


1 3 

A m(h) = hbo+a^h^F 




h 2 +vi{m— \)h*+V 2 {m—\) 2 h?+0{tU ), 


uniformly for h £ (0, ho) and m £ Nf )]C . 

In particular, the splitting between two consecutive eigenvalues satisfies 


\ m+1 {h)-\ m {h) = dh 2 + o {hi). 


Proof. If the support of k * is small enough, the hypothesis k*{h^,Z) = 

3 

<D((h + |Z|)2) implies that, when h is small enough, 

2 1 

(1 + z?)/C^ ^ K-n + -k*(hz , ICji) ^ (1 — r/)ICfi , 

for some small 77 > 0. Therefore, since the eigenvalues of JCn are (2 m — 1 )h, 
m £ N*, the variational principle implies that the number of eigenvalues of 
K h +lk*{hh,K h ) below a threshold belongs to f^ +1), 

Taking Cn = ^{c — 1 )o l t 2 h 1 / 2 + ^h, and applying the theorem, we obtain 
the estimate for the cardinal of N^ c . The corresponding eigenvalues of Cn,A 
are of the form 


1 3 C 

A m {h) = hb 0 + aihi - {Ajh 2 + h 


0{m-l) + k*{hi,2m-l)} +0(h°°), 
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with (2m — 1 )h ^ Therefore there exists a constant C > 0, independent 
of h, such that all m G IM/^ satisfy the inequality (2m — l)h ^ Ch 1//2 . Writing 


k*(h^,Z) = c 0 h 3/2 + vih l/2 {Z/2) + cih 2 + v 2 {Z/2) 2 + v 3 hZ 

+ h l/2 0(h + \Z\) 2 + 0{Z 3 ), 


we see that, for m G N^ iC , 

m 2 J 

which gives the result. □ 


k*(f 12 , (2m — l)/i) = v\h 3 / 2 ^ 



Remark 1.6. An upper bound of A m {K) for fixed ^-independent m with re- 

9 

mainder in 0{h 4 ) was obtained in [13] through a quasimodes construction 
involving powers of h 4. To the authors’ knowledge, Corollary 1.5 gives the 
most accurate description of magnetic eigenvalues in three dimensions, in 
such a large spectral window. Note also that the non-degeneracy assump¬ 
tion on the norm of B is not purely technical. Indeed, at the quantum level, 
it appears through microlocal reductions matching with the splitting of the 
Hamiltonian dynamics into three scales: the cyclotron motion around field 
lines, the center-guide oscillation along the field lines, and the oscillation 
within the space of field lines. 


1.6. Organization of the paper. The paper is organized as follows. In 
Section 2, we state our main results. Section 3 is devoted to the investigation 
of the first normal form (see Theorem 2.1 and Corollary 2.4). In Section 4 we 
analyze the second normal form (see Theorems 2.8 and 2.11 and Corollaries 
2.9 and 2.13). Section 5 is devoted to the third normal form (see Theorem 
2.15 and Corollary 2.16). 


2. Statements of the main results 

We recall (see [9, Chapter 7]) that a function m : —> [0, 00 ) is an order 

function if there exist constants Nq, Cq > 0 such that 

m(X) ^ C 0 (X - Y) N °m(Y) 

for any X,Y G M rf . The symbol class S(m ) is the space of smooth h- 
dependent functions —> C such that 

Va G N d , |5“a^(x)| ^ C a m(x), V/i G (0,1]. 

Throughout this paper, we assume that the components of the vector 
potential A belong to a symbol class S(m). Note that this implies that 
B G S(m), and conversely, if B G S(m), then there exist a potential A 
and another order function ml such that A G S(m'). Moreover, the magnetic 
Hamiltonian H(x,£) = ||£ — A(x)|p belongs to S(m") for an order function 
m" on M 6 . 
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We will work with the Weyl quantization; for a classical symbol an = 
a(x,£;h) E S(m), it is defined as: 

( 2 . 1 ) 

Op™ aiffx) = 1 [ e l{x ~ y ^ /h a ip(y)dyd£, \/if E S(R d ). 

(2TTh) d J R 2d V 2 ) 

The Weyl quantization of H is the magnetic Laplacian Cn ,a = (—iKV — A) 2 . 

2 . 1 . Normal forms and spectral reductions. Let us introduce our first 
Birkhoff normal form A//,,. 

Theorem 2.1. //B(0) 0, there exists a neighborhood of { 0, A(0)) endowed 

with symplectic coordinates (x±, £i, X 2 , £ 2 , % 3 , £ 3 ) in which £ = {aq = £1 = 
£3 = 0} and (0, A(0)) has coordinates 0 E M 6 , and there exist an associated 
unitary Fourier integral operator Un and a smooth function, compactly sup¬ 
ported with respect to Z and £3, f*(h,Z,x 2, £2^3) £3) whose Taylor series 
with respect to Z,^,h is 

k ^3 2(+2m+l3=k 

such that 

( 2 . 2 ) U* h C h ,AU h =M h + lZ h , 

with 

M h = h 2 Dl 3 + In Op^ b + Op™ f*{h,l h ,x 2, £2,^3, £3), 

and where 

(a) we have In = H 2 D 2 1 + x\, 

(b) the operator Op™ f*(h,ln,x 2, £2, £3, £3) has to be understood as the Weyl 
quantization of an operator valued symbol, 

(c) the remainder IZn is a pseudo-differential operator such that, in a neigh¬ 
borhood of the origin, the Taylor series of its symbol with respect to 
(a=i,£i,£3 ,K) is 0 . 

Remark 2 . 2 . In Theorem 2 . 1 , the direction of B considered as a vector field 
on £ is and the function b E C°°(M 6 ) stands for b o o it where it : 

M 6 —> £ : tt(x±, £1, X2, £2, %3, £3) = (0,0, X2, £2, X3, 0). In addition, note that 
the support of /* in Z and £3 may be chosen as small as we want. 

Remark 2 . 3 . In the context of Weyl’s asymptotics, a close version of this 
theorem appears in [ 19 , Chapter 6 ]. 

In order to investigate the spectrum of Cn,A near the low lying energies, 
we introduce the following pseudo-differential operator 

a4 1] = h 2 Dl 3 + h Op£ b + 0 P ;r f(h, h, ® 2 , 6 , * 3 , 6 ), 
obtained by replacing X/ ( by h. 
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Corollary 2.4. We introduce 

(2-3) JVj = CW (iVjj) , 

with 

N h = €3 + I hb(x 2 , 6 ^ 3 ) + f*’\h,T n ,X 2 ,&,X 3 ,f 3 ) 
and where b is a smooth extension of b away from D{ 0,e) such that (1.11) 
still holds and where /*’** = x( x 2 ,&, x 3 )/*, with x is a smooth cutoff function 
being 1 in a neighborhood of D(0,e). We also define the operator attached to 
the first eigenvalue ofXji 

(2.4) A/J 1] ’ # = Op£ (l\f ] ’“) , 

where N™ = + ^k{x 2 ,&,x 3 ) + f* A {K K % 2 , &,x 3 , &)■ 

If e and the support of f* are small enough, then we have 

(a) The spectra of a and below coincide modulo 0(h°°). 

r-, 1 u 

(b) For all c E (0, min(36o, fio)), the spectra of £/,,a and A/)( below ch 
coincide modulo 0 (K°°). 

Let us now state our results concerning the normal form of (or A/l*^') 
under the following assumption. 

Notation 2.5. If / = /(z) is a differentiable function, we denote by T z /(-) 
its tangent map at the point z. Moreover, if / is twice differentiable, the 
second derivative of / is denoted by T %/(■, •). 

Assumption 2.6. We assume that Tq 6 (B( 0 ), B(0)) > 0. 

Remark 2.7. If the function b admits a unique and positive minimum at 0 
and that it is non degenerate, then Assumption 2.6 is satisfied. 

Under Assumption 2.6, we have d 3 b(0, 0,0) = 0 and, in the coordinates 
( x 2 , £, 2 , X 3 ) given in Theorem 2.1, 

(2.5) 6(0,0,0) > 0. 

It follows from (2.5) and the implicit function theorem that, for small x- 2 , 
there exists a smooth function (x 2 , £ 2 ) >->• s(x 2 , £ 2 ), s( 0 , 0 ) = 0 , such that 

( 2 . 6 ) d 3 b(x2,f2is(x2,&)) = 0 . 

The point s(x 2 , £ 2 ) is the unique (in a neighborhood of (0, 0,0)) minimum of 
x 3 (-)• b(x 2 ,^ 2 ,x 3 ). We define 

^(* 2 , 6 ) := (5dlb(x2,&,s(x 2 ,&))) 1/4 - 

Theorem 2.8. Under Assumption 2.6, there exists a neighborhood Vo of 0 
and a Fourier integral operator Vn which is microlocally unitary near Vo and 
such that 

V h *M [4] V n =:^ ] = Op£ (4 1] ), 


12 B. HELFFER, Y. KORDYUKOV, N. RAYMOND, AND S. VU NGOC 

where = u 2 (x 2 ,f 2 ) (£3 + hxf) + hb(x 2 ,£ 2 ,s(x 2,6)) + Ln an d Ln is a 
semiclassical symbol such that r h = (D(hx\) + 0 (fU ; 2 ) + 0(£f) + 0 {h 2 ). 

Corollary 2.9. Let us introduce 

JSLf A = Op^ (4 1] ’“) » 

where = ^ 2 (x 2 ,6) (£3 + hxf) + hb(x 2 , £ 2 , s(x 2 , £ 2 )) + with r\ = 

x(®2, £2, £3, &)Zh> and where 0 denotes a smooth and constant (with a posi¬ 
tive constant) extension of the function v. 

There exists a constant c > 0 such that, for any cut-off function x equal to 
1 on D{ 0, e) with support in D{ 0, 2e), we have: 

(a) The spectra of A/j- 1 ^ andj\f\below (bo + ce 2 )h coincide modulo 0(h°°). 

(b) For all c E (0, min(3&0i &o + ce 2 )), the spectra of Cn ,a and below ch 
coincide modulo 0{h°°). 

Notation 2.10 (Change of semiclassical parameter). We let h = h? and, if 
An is a semiclassical symbol on T*M 2 , admitting a semiclassical expansion in 
, we write 

A h := Op^ A h = Opf A h =: ?l h , 

with 

&h{x 2 ,£ 2 ,x 3 ,£ 3 ) = A h 2 (x 2 ,h£ 2 ,x 3 ,h£ 3 ). 

Thus, An and 21^ represent the same operator when h = h ?, but the former 
is viewed as an ^-quantization of the symbol An, while the latter is an h- 
pseudo-differential operator with symbol A/,.. Notice that, if An belongs to 
some class S(m), then A n E S(m) as well. This is of course not true the other 
way around. 

Theorem 2.11. Under Assumption 2.6, there exist a unitary operator Wh 
and a smooth function g*{h , Z,x 2 ,f 2 ), with compact support as small as we 
want with respect to Z and with compact support in (x 2 ,£ 2 ), whose Taylor 
series with respect to Z, h is 

J2 C mAx2 ,&)Z m h e , 

2 m+ 21>3 

such that 

W* h fi U W h =-.m h = 0 p)"(M fc ), 

with 

Mh = h 2 b(x 2 , h£ 2 , s(x 2 , hf 2 )) + h 2 J h Op™ n 2 (x 2 , h£ 2 ) + h 2 g*(h, J h , x 2 , h £ 2 ) 

+ h 2 R h + h°°S(l). 

where 

fii u r-i I u 

(a) the operator y\ l h ' 9 is AT) (but written in the h-quantization), 

(b) we have let Jn = Opf (j) 2 + x^j, 
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(c) the function R h satisfies R h (x 2 , h£ 2 , x 3 , £ 3 ) = 0((x 3 ,£ 3 )°°). 

Remark 2.12. Note that the support of g * with respect to Z may be chosen 
as small as we want. Note also that we have used 91^'^ instead of 91^: Since 
Wh is exactly unitary, we get a direct comparison of the spectra. 

Corollary 2.13. We introduce 

K = CW (m») , 

with 

= h 2 b(x 2 ,h£ 2 ,s(x 2 ,h £ 2 )) + h 2 J h u 2 (x 2 ,h£ 2 ) + h 2 g*(h,J h ,x 2 ,h£ 2 ). 

We also define 

< 1] ’“ = Ovl (mW’“) , 

with 

M™ = h 2 b(x 2 , h£ 2 , s(x 2 , h£ 2 )) + h 3 n 2 (x 2 , hf 2 ) + h 2 g*(h , h, x 2 , h£ 2 ). 

If e and the support of g* are small enough, we have 

(a) For all g > 0, the spectra o/9t^ and 9Jl| below boh 2 + 0(h 2+ri ) coincide 
modulo 0(h°°). 

(b) Fore 6 (0,3), the spectra o/9Jt^ and below boh 2 + ca^h? coincide 

modulo 0(h°°). 

(c) If c € (0,3), the spectra of a and M.\^' = 911 below boh + ccr^h? 
coincide modulo 0{h°°). 

r-i I n 

Finally, we can perform a last Birkhoff normal form for the operator My h * 
as soon as (X 2 , £ 2 ) ^ H x 2 i^ 2 , s(x 2 ,f 2 )) admits a unique and non degenerate 
minimum at (0, 0). Under this additional assumption, b admits a unique and 
non degenerate minimum at (0,0,0). 

Therefore we will use the following stronger assumption. 

Assumption 2.14. The function b admits a unique and positive minimum 
at 0 and it is non degenerate. 

Theorem 2.15. Under Assumption 2. If, there exist a unitary h-Fourier 
Integral Operator Q 1 whose phase admits an expansion in powers of h 2 
such that 

where 

(a) J-fr is defined in Theorem 1.3, 

(b) the remainder is in the form Q% = Op^’ ( Gn), with Gn = hO(\z 2 \°°). 

Corollary 2.16. If e and the support of k* are small enough, we have 

(a) For all g G (0, ^), the spectra of and Fh below boh + 0(h 1+v ) 

coincide modulo 
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(b) For all c E (0, 3), the spectra of Cn,A arid below b^h + ca^h? coincide 
modulo 0(h°°). 

Remark 2.17. Since the spectral analysis of J-ji is straightforward, Item (b) 
of Corollary 2.16 implies Theorem 1.3. 

The next sections are devoted to the proofs of our main results. 


Corollary 2 . 4 b 




3. First Birkhoff normal form 


We assume that B(0) 7 ^ 0 so that in some neighborhood 0 of 0 the mag¬ 
netic field does not vanish. Up to a rotation in M 3 (extended to a symplectic 
transformation in M 6 ) we may assume that B(0) = ||B(0)||e 3 . In this neigh¬ 
borhood, we may defined the unit vector: 


(3.1) 



and find vectors c and d depending smoothly on q such that (b, c, d) is a 
direct orthonormal basis. 


3.1. Symplectic coordinates. 

3.1.1. Straightening the magnetic vector field. We consider the form do and 
we would like to find a diffeomorphism, in a neighborhood of 0 , x such that 
x(q) = q and x*(da) = dgi A dq 2 ■ First, this is easy to find a local diffeomor¬ 
phism tp such that 

d 3 <p(q) = b (ip(q)) 

and ip(qi,q2,0) = (gi, <? 2 ,0). This is just the standard straigthening-out 
lemma for the non-vanishing vector field b. 

The vector e 3 is in the kernel of ip* (da), which implies that we have 
yA(da) = f(q)dq\ A dq 2 , for some smooth function /. 
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But since the form p*(da) is closed, / does not depend on q 3 . This is 
then easy to find another diffeomorphism x/j, corresponding to the change of 
variables 

q = ^(q) = £ 2 ),^ 2 ( 91 , £ 2 ), &), 

such that 

ip* (if * (da)) = dqi A dq 2 • 

We let x = P 0 V' and we notice that 

(3.2) X*(da) = dqi A d q 2 d 3 x(q) = Hx(q)) , 

Remark 3.1. It follows from (3.2) and (1.4) that detTx = ||B|| -1 . 

3.1.2. Symplectic coordinates. Let us consider the new parametrization of E 
given by 

i : n —y E 

q (x(q),Ai(x(q )), A 2 (x(q)),A 3 ( X (q ))), 
which gives a basis (fi, f 2 , fa) of T E : 

fj = ( T x(ej), T A o Tx(ej)), j = 1, 2,3 . 

Using (1.5), and the fact that f 3 is in the kernel of da, we find wo(f?, fa) = 0, 
j = 1,2. Finally, u 0 (fi,f 2 ) = da^Tyex, Tye 2 ) = x*(da)(e 4 , e 2 ) = 1. 

The following vectors of M 3 x M 3 form a basis of the symplectic orthogonal 

°f T d?)E: 

(3.3) f 4 = llBir 1 / 2 ^, ( f T x (flA)c), f 5 = ||B||- 1 / 2 (d, (*T xW) A)d), 
so that 

Wo(f4, f5) = —1. 

We let fg = (0, b) + pifi + /? 2 f 2 where p\ and p 2 are determined so that 
cc 0 (fj, fg) = 0 for j = 1,2. We notice that cuo(fj; fg) = 0 for j = 4,5 and 
Wo(f3j fb) = —1- 

3.1.3. Diagonalizing the Hessian. We recall that 

H(q,p) = \\p- A(q)\\ 2 

so that, at a critical point p = A (q), the Hessian is 

T 2 H((U u Vi), (U 2 , V 2 )) = 2(Ui - T q A(U{), V 2 - T q A(U 2 )). 

Let us notice that 

T 2 H( f 4 ,f 5 ) = 21|B || —1 (B x c,B x d) = 0, 

T 2 iL(f 4 , fg) = 2(B x c, b) = 0, 

T 2 H(i 3 , fg) = 2(B x d, b) = 0. 

The Hessian is diagonal in the basis (f 4 , fg, fg). Moreover we have 
r 2 iL(f 4 ,f 4 ) = d 2 iL(f 5 ,f 5 ) = 2||B|| _ 1 1|B x c || 2 = 2||B||~ 1 ||B x d|| 2 = 2||B||. 
Finally we have: 


T 2 H( f 6 ,f 6 ) = 2. 
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Now we consider the local diffeomorphism: 

(x,£) H- i{x 2 ,&,X 3 ) + Xif 4 (x 2 ,£ 2 ,X 3 ) + ilh{.X 2 ^ 2 ,X^) + £sfe(x 2 , £ 2 , X 3 ). 

The Jacobian of this map is a symplectic matrix on X. We may apply the 
Moser-Weinstein argument (see [29]) to make this map locally symplectic 
near X modulo a change of variable which is tangent to the identity. 

Near X, in these new coordinates, the Hamiltonian H admits the expan¬ 
sion 

(3.4) H = H° + 0(M 3 + |£i | 3 + |£ 3 | 3 ), 

where H denotes H in the coordinates (xi, X 2 , X3, £1, £ 2 , £ 3 ), and with 

(3.5) H° = £f + &(x 2 ,£ 2 ,x 3 )(x? + £?), b = ||B(x 2 ,&,x 3 )||. 

3.2. Semiclassical BirkhofF normal form. 

3.2.1. Birkhoff procedure in formal series. Let us consider the space £ of 
formal power series in (xi,£i,£ 3 , h) with coefficients smoothly depending on 
x = (x 2 ,6,x 3 ): 

We endow £ with the semiclassical Moyal product (with respect to all vari¬ 
ables (xi, x 2 , x 3 , £ 1 , £ 2 , £ 3 )) denoted by * and the commutator of two series K\ 
and K 2 is defined as 

[«i, k 2 ] = * k 2 - H 2 * «i . 

The degree of = zf^h 1 is a\ + a 2 + P + 2-f = |a| + f3 + 21. Vn 

denotes the space of monomials of degree N. On is the space of formal series 
with valuation at least N. For any r, 7 G £, we denote ad T 7 = [r, 7]. 

Proposition 3 . 2 . Given 7 G O3, there exist formal power series r, n G O3 
such that 

e ift- 1 adT( F o + 7 ) = h o + Ki 

with [k, \z \| 2 ] = 0 . 

Proof. Let N ^ 1. Assume that we have, for tn G O3, 

e^ 1 ad ^v (H° + 7 ) = h° + K 3 + --- + K n+ 1 + R n+2 + O n+ 3 , 

with Ki G Vi, [K i: \zi\ 2 } = 0 and Rn+2 € V N +2 ■ 

Let t' G Vn+ 2 - Then we have 

e^' 1 ad (H° + 7 ) = H° + K 3 + • • • + K n+1 + K n +2 + O n+3 , 

with Kn+2 G Vn+2 such that 

Kn +2 = Rn +2 + ih 1 ad T ' iL° + Otv+ 3 - 
Lemma 3 . 3 . For t' G Vn+2; ?xe Ziaxe 

ih ~ l ad T ' H° = ih~ 1 b ad r / |<;i | 2 + 0 at +3 . 
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To prove this lemma, we observe that 

ih -1 ad T / H° = ihr 1 ad T / + ihT 1 ad T /(&(x)|zi| 2 ). 


Let us write 

T = a Q>( j t i(x)z^§n l . 

\a\+/3+2l=N+2 

Then, for the first term, we have 
1 c2 _r > t:2^ 


ih ad r / £ 3 ={t , £ 3 } 

= -26§^ 

dx 3 

= -* E 


\a\+/3+2£=N+2 


d a a,P ,i (~\„a C P +1 
dx 3 


(x)zur 


We also have 

ih~ 1 ( ad T / 6(x)) ={r / , b} + K 2 On 

dr' db dr' db dr’ db 

<9£3 <9 x '3 + <9£ 2 <9^2 5^2 <96 + ^ +1 

c)b 

E + Oiv +1 € o N+1 • 

\a\+/3+2£=N+2 

Therefore, for the second term, we get 

ih -1 ad r /(6(x)|zi| 2 ) =ih~ 1 ( ad T / 6(5;))|^i| 2 + ih~ 1 b(x) ad T / |zi| 2 

=ih~ 1 b(x) ad T / l^i | 2 + On +3 , 

that completes the proof of the lemma. 

By the lemma, we obtain that 

Kn+2 = Kn+2 + ba.d T i \zi | 2 , 


that we rewrite as 

Rn+2 = Kn+2 + ih ^adi^i 

Since b(x) / 0, we deduce the existence of r l and I\n +2 such that Kn +2 
commutes with |zi| 2 . □ 


2 t = K n +2 + 6{|^i| 2 , t’}. 


3.2.2. Quantizing the formal procedure. Let us now prove Theorem 2.1. Us¬ 
ing (3.4) and applying the Egorov theorem (see [26, 31] or Theorem A.2), we 
can find a unitary Fourier Integral Operator Un, and such that 

U; t C h .AU h = C 0 h + Op %(H°) + Op %(r n ), 


where the Taylor series (with respect to x\, £1, £3, h) of r^ satisfies rj = 
7 E O 3 and Cq is the value at the origin of the sub-principal symbol of 
One can choose U^ such that the subprincipal symbol is preserved 
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by conjugation 1 , which implies Co = 0. Applying Proposition 3.2, we obtain 
r and n in O 3 such that 

e ih-^d T ^ H o + 7 ) = H° + K, 

with [k, \zi\ 2 ] = 0. 

We can introduce a smooth symbol an with compact support such that 
we have af = r in a neighborhood of the origin. By Proposition 3.2 and 
Theorem A.4, we obtain that the operator 

e ih ~ 10p n( a *-\Op%(H 0 ) + Op %(r h ))e- ih ~ l0p * M 

is a pseudodifferential operator such that the formal Taylor series of its sym¬ 
bol is H° + k. In this application of Theorem A.4, we have used the filtration 
Oj defined in Section 3.2.1. Since k commutes with |zi| 2 , we can write it as 
a formal series in \z\\ 2 : 

K = ^2 ^ c i,m{x2^2,X3)\z l \ 2m ^. 

k^3 2 l+ 2 m+P=k 

This formal series can be reordered by using monomials {\z\\ 2 )* m : 
k = J2 hic lm(x 2 ,&,X 3 ){\z 1 \ 2 )* m 

k^ 3 2i+2m+f3=k 

Thanks to the Borel lemma, we may find a smooth function, with a compact 
support as small as we want with respect to h, I and £3, f*(h, /, X2, £ 2 , £3, £3) 
such that its Taylor series with respect to h , /, £3 is 

k^3 2i+2m+/3=k 

This achieves the proof of Theorem 2.1. 

3.3. Spectral reduction to the first normal form. This section is de¬ 
voted to the proof of Corollary 2.4. 

3.3.1. Numbers of eigenvalues. 

Lemma 3.4. Under Assumption 1.2, there exists ho > 0 and £0 > 0 such 
that for all h £ ( 0 , ho), 

infSess(A^) ^ (/3o + e 0 )h. 

Proof. By using the assumption we may consider a smooth function x with 
compact support and £o > 0 such that 

£3 + b( x 2 , £ 2 , x 3 ) + x( x 2 , x 3l £ 2 , £ 3 ) > A) + 2 eo- 


1 This is sometimes called the Improved Egorov Theorem. It was first discovered by We¬ 

instein in [30], in the homogeneous setting. For the semiclassical case, see for instance [18, 
Appendix A]. 
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Then, given r/ E (0,1) and estimating the second term in (2.3) by using that 
the support of /* is chosen small enough and the semiclassical Calderon- 
Vaillancourt theorem, we notice that, for h small enough, 

(3.6) Ml ^ (1 - J 7 )Op)? (£§ + \z 1 \ 2 b{x 2 ,C 2 ,x 3 )) . 

Since the essential spectrum is invariant by (relatively) compact perturba¬ 
tions, we have 

Sess (aC| + (1 - r/)/i Op)" x(x 2 , £3,6, £ 3 )) = Sess (.A/j) . 

Hence 

infSess (a/«) ^ infs (a(\ + (1 - rfjhOp™ x(x 2 , x 3 , £2, £3)) • 

In order to bound the r.h.s. from below, we write 

A/jj + (l-#Op^x(a;2,a:3,£2,£3) 

> 0-~V)OPh (£3 + |^l| 2 6(x 2 ,£2,.'C 3 )) + (1 -I/)^Op^x(^2,a;3,£2,£3) 

> ^(1 -1\) Op“ (£3 + b(x 2 , £ 2 , x 3 ) + X( x 2 , X 3 , £ 2 , £ 3 )) 

^ ft(l - r/)(/ 3 0 + 2£0 - CTt), 

where we have used the semiclassical Garding inequality. Taking rj and then 
h small enough, this concludes the proof. □ 

By using the Hilbertian decomposition given by the Hermite functions 
{&k,h)k^i associated with Z/,, we notice that 



k^i 


where 

( 3 . 7 ) M^ = h 2 D 2 X3 + (2k-l)hOp^b+Op^ f^(h,(2k-l)h,x 2 M2,x 3 ,^), 
acting on L 2 (M 2 ). 

Lemma 3.5. For all ?y E (0,1), there exist C > 0 and ho > 0 such that for 
all k ^ 1 and h E (0, ho), we have Si(M^’^) ^ (1 — 2g)bo(2k — 1 )h. 

Proof. Applying ( 3 . 6 ) to V>(xi, x 2 , x 3 ) = 99 ( 0 : 2 , x 3 )^k,h( x i)> we infer that 

(AT^’V,^) > (2A: - l)ft(l - ??) (Op™ (%,</>)• 

With the Garding inequality, we get 

(OP 'n(b)<p,<p) ^ ( b 0 - Ch)\\<p\\ 2 , 

and the conclusion follows by the min-max principle. □ 

We immediately deduce the following proposition. 

Proposition 3.6. We have the following descriptions of the low lying spec¬ 
trum of Mf). 
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(a) There exist ho > 0 and K G N such that, for h E (0, ho), the spectrum of 
A lying below /3 q h is contained in the union UfcLi s p (j^h U ) ■ 

(b) If c E (0, min(36o, (do)), then there exists ho > 0 such that for all h E 
(0, ho) the eigenvalues of Aft lying below ch coincide with the eigenvalues 
of below ch. 

We deduce the following proposition. 

Corollary 3.7. Under Assumption 1.10, we have 

N (£ n ,A,Poh) = o(h~ 3 / 2 ), N (a f*,Poh) = o(h~ 2 ). 


Proof. To get the first estimate, we use the Lieb-Thirring inequalities (which 
provide an upper bound of the number of eigenvalues in dimension three) 
and the diamagnetic inequality (see [25] and Proposition 1.8). To get the 
second estimate, we use the first point in Proposition 3.6. Moreover, given 
r] E (0,1), by using h E (0,1) we infer 

(A ^ (1 - r])h(0 p™ (C| + b(x 2 ,f 2 ,x 3 )) 

Note that the last inequality is very rough. By the min-max principle, we 
deduce that 


N A 


[kU Potty < N (Op£ (C| + b(x 2 , C 2 , ®3)) ,(!-»/) Vo) • 

Then, we conclude by using the Weyl asymptotics and our confinement as¬ 
sumption: 

N (Op™ (C| + b(x 2 ,&,x 3 )) , (1 - h)~ 1 Po) = 0(h~ 2 ). 

□ 

Since A/]? commutes with X/,, we also deduce the following corollary. 

Corollary 3 . 8 . For any eigenvalue A of A/]? such that A ^ PqK we may 
consider an orthonormal eigenbasis of the space ker ^A/"| — A^ formed with 
functions in the form ft(xi)<p%(x 2 , £3) with k E {1 Moreover we 

have 1(-oo j/ 3oM (A/)?) = 0(fh~ 2 ) and each eigenfunction associated with A ^ 
Poh is a linear combination of at most 0 (h~ 2 ) such tensor products. 

3.3.2. Microlocalization estimates. The following proposition follows from 
the same lines as in dimension two (see [16, Theorem 2.1]). 

Proposition 3 . 9 . Under Assumptions 1.1 and 1.2, for any e > 0, there 
exist C{e) >0 and ho(e) > 0 such that for any eigenpair (A, if) of Cn, a with 
A ^ Pq h we have for h E (0, /10(e)): 

f e 2 il ~ e) ^ q)/hh \if \ 2 dq ^ C(e)exp(e?T*] 


Q/i,A(e (1 ' e)0(9)// ^V < C(e)exp(e7T5 
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where <f is the distance to the bounded set {||B(g)|| ^ (3 0} for the Agmon 
metric (||(B(g)|| — f3o)+g, with g the standard metric. 

Proposition 3.10. Under Assumptions 1.1 and 1.2, we consider 0 < bo < 
/3o < b\ and there exist C > 0 and ho > 0 such that for any eigenpair (A, if) 
of Cn,A with A ^ /3oh we have for h 6 (0, ho) and 6 6 (0, ^): 

= Xo (h~ 2S Cn, a) XI( 9)^ + o(h°°)Ul 

where xo a cutoff function compactly supported in the ball of center 0 and 
radius 1 and where xi is a compactly supported smooth cutoff function being 
1 in an open neighborhood o/{||B(g)|| ^ /3q}. 

Let us now investigate the microlocalization of the eigenfunctions of Aft. 

Proposition 3.11. Let x be a smooth cutoff function being 0 on {6 ^ /3o} 
and 1 on the set {b ^ (3o + e}. If A is an eigenvalue o/A/"| such that A ^ /3oh 
and if if is an associated eigenfunction, then we have 

Op™ (x(*2,6,*3 ))if = o(h°°)\\if\\. 


Proof. Due to Corollary 3.8, it is sufficient to prove the estimate for a func¬ 
tion in the form if(xi,X 2 ,x 3 ) = ek,h(xi)if(x 2 , xf) where k lies in {1,, I\ } 
and we have 

Aflif = \if , or equivalently = A ip, 

where we recall (3.7). Then, we write 


A/fOp)T (X) <P = A OpJ" (X) + H kU , Op£ (X) 


<p 


and it follows that 


(3.8) 


(a f u <W (x) V, Ovt (x) f) = All Op}» (x) vtf 


+ 



Rough pseudo-differential estimates imply that there exist C > 0, ho > 0 
such that for all h £ (0, ho), 


(3.9) 


Opj? (x) <p, Opf (x) pj < Ch 2 ||Opjf (x) p\\ 2 +ch ||o P £ (x) <p\ 
+ Ch(Op™ (d 3X ) <p, Op% (&) o P j? (x) ip). 


Combining (3.9) and (3.8), we get 

(3.10) 11 Op™ (£3) Opjf (x) V\\ < Ch 5 || Op£ (x) ip\\ , 


where x is a smooth cutoff function living on a slightly larger support than 
X- By using (3.10), we can improve the commutator estimate 


A/f J ,Op£ (x)l <P,Op% (x) P 


€ 


Ch ~2 


I°p^ (x) A 
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We infer that, there exist C > 0, Hq > 0 such that for h E (0, Ho), 

(wf 1,# °P n (x) <P, °Ph (x) ¥>) < W\ Op n (x) <p\\ 2 + cn ^ ||Op“ (x) pf ■ 

By using the semiclassical Garding inequality and the support of y, we get 

(4 kU °Pfr (x) <P, Opjf (x) <p) > (A) + 11 Op™ (x) v\\ 2 

and we deduce 

IIOpr (x) <p\\ 2 < cnI ||Opjr (x) <^|| 2 • 

The conclusion follows by a standard iteration argument. □ 


The following proposition is concerned by the microlocalization with re¬ 
spect to £3. 

Proposition 3.12. Let xo be a smooth cutoff function being 0 in a neigh¬ 
borhood of 0 and let 5 E (0, |) • If X is an eigenvalue of such that A fioh 
and if if is an associated eigenfunction, then we have 

Opf (xo (h- 5 ^if = o(h°°)\\if\\. 

Proof. We write again if(xi,x 2 ,x 3 ) = ek,h(xi)<p(x 2 , x 3 ) with k E {1,..., K} 

and we have A4^’V = We use again the formula (3.8) with xo (^ _<5 £ 3 ). 

We get the commutator estimate 


A/f J ,Op£ (xo (h 5 &)) <P, Op™ (xo (h ^3)) <P 

^ ||o P ^ (x 0 (ft- 5 £ 3 ))^ 


We have 


Op™ ({h 5 6) 2 X'o *&)) = Op£_, (?lxo (6)) , 

so that, with the Garding inequality, 


Op£ ((/t^ 3 ) 2 xo (/t*£ 3 )) <p,<p) > (1 - Ch l ~ 5 


We infer 


[h 25 { 1 - Ch 1 - 5 ) - Poh) 

Op™ (xo (h 5 £ 3 )) P 

2 


< Chi~ 5 

Op H (x 0 ( h ~ S &) ) <P 


□ 


Using Op% f*(h,l h ,x 2 ,£2,X3,£ 3 ) = Op% f(h, \z 1 \ 2 ,x 2 ,&,X3,&), we de¬ 
duce the following in the same way. 

Proposition 3.13. Let xi be a smooth cutoff function being 0 in a neigh¬ 
borhood of 0 and let 5 E (0, . If A is an eigenvalue of Mf such that A ^ /3q h 

and if if is an associated eigenfunction, then we have 

Op™ (xi (/T*(xi,fi))) if = 0(h°°)\\if\\. 
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Proposition 3.14. The spectra of Cn,A and J\f^ below /3oh coincide modulo 

Proof. We refer to [25, Section 4.3] which contains similar arguments. □ 

This proposition provides the point (a) in Corollary 2.4. With Proposi¬ 
tion 3.6, we deduce the point (b). 

4. Second Birkhoff normal form 

4.1. Birkhoff analysis of the first level. This section is devoted to the 
proofs of Theorems 2.8 and 2.11. 

The goal now is to normalize a li-pseudo-differential operator A/"^ on M 2 
whose Weyl symbol has the form 

+ hb(x 2 ,&,x 3 ) +r h (x 2 ,&,x 3 ,£ 3 ), 

where rn is a classical symbol with the following asymptotic expansion: 

rn = to + hri + h 2 r 2 H- 

(in the symbol class topology), where each ri has a formal expansion in £3 of 
the form 

(4.1) r e (x 2 ,&,x 3 ,£ 3 )~ c iA x 2^2, x 3)^ 3 - 

2 £+/ 3^3 

The leading terms of are: 

(4.2) ivf 1 =^+hb(x 2 ^ 2 ,x 3 )+c lA (x 2 ,^,x 3 M 3 +0(h^i)+0(e 3 )+0(h 2 ). 

4.1.1. First normalization of the symbol. We consider the following local 
change of variables <p(x 2 , £ 2 , x 3 , £3) = (x 2 , £ 2 , x 3 , £3): 

' x 2 := x 2 +f 3 d 2 s(x 2 ,f 2 ), 

^ ^ 6 := 6 + 6^is(®2,6) , 

£3 ■= x 3 - s(x 2 ,f 2 ), 

„ £3 := £3 • 

It is easy to check that the differential of (p is invertible as soon as £3 is small 
enough. Moreover, we have 

(p*Uo — w 0 = O (I ^31 ). 

By the Darboux-Weinstein theorem (see for instance [25, Lemma 2.4]), there 
exists a local diffeomorphism if such that 

(4.4) if = Id + £}(£!) and cao = wo- 

Using the improved Egorov theorem, one can find a unitary Fourier Inte¬ 
gral Operator V% such that the Weyl symbol of V^Af^Vn is Nn := o o 
if + 0(h 2 ). From (4.4), and (4.3), we see that f% := r^o (p o if is still of the 
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form ( 4 . 1 ), with modified coefficients Cf.g. Thus, using the new variables and 
a Taylor expansion in £3, we get 

N h = il + hb(x 2 + 0 (t 3 ),£ 2 + 0 (&),x 3 + s(x 2 + 0 (i 3 ),i 2 + 0 (t 3 )) + 0 (il)) 

+ 0(£f) + + 0(h 2 ) 


and thus 

( 4 . 5 ) N h = £| + hb(x 2 , £2, x 3 + s(x 2 , £2)) + hi 3 g(x 2 , £2, £3) 

+ I) + n* + cK£l) + 


for some smooth function g{x 2, £2, £3). 

Therefore A), has the following form: 

= £|+^6(x 2 , £2, x 3 +s(x 2 , £2 ))+ci,i(x 2 , £2, x 3 )/i£3+C>(^£|)+C>(£|)+C>(^ 2 ). 

4 . 1 . 2 . W/iere t/ie second harmonic oscillator appears. We now drop all the 
hats off the variables. We use a Taylor expansion with respect to x 3 , which, 
in view of (2.6), yields: 


fr( x 2, £2, x 3 + s(x 2 , £ 2 )) = 6(®2, £2, s(®2, £2)) + y 6(x 2 , £2, a(x 2 , £2)) + 0 (x|). 
We let: 

(4.6) v = ( 23 | 6 (x 2 ,£ 2 ,s(x 2 ,£ 2 ))) 1/4 and 7 = In ;y. 

We introduce the change of coordinates (x 2 , X 3 , £2, £3) = C(x 2 , £3, £2, £3) de¬ 
fined by: 

= i/x 3 , 

,, £3 = i /_1 £3, 

^ ^ x 2 = x 2 + §^x 3 £, 3 , 

k £2 = £2 - £®3£3 , 

for which one can check that C*ujq — ujq = 0(x 3 £ 3 ) = 0(£ 3 ). As before, we 
can make this local diffeomorphism symplectic by the Darboux-Weinstein 
theorem, which modifies (4.7) by C?(£|). In the new variables (which we call 
(X 2 , X 3 , £ 2 , £3) again), the symbol Nn has the form: 

N h = V 2 (X2,&) (£3 +hx\) + hb(x 2 , £2, s(x2, £2)) +Cl i l(x2,£2,X 3 )^£ 3 

+ o(hxl) + <D(h$) + o (£ 3 3 ) + 0(h 2 ), 


for some smooth function ci.i(x 2 , £2, x 3 ). 
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4.1.3. Normalizing the remainder. The next step is to get rid of the term 
ci,i(x 2 ,&,x 3 )H £3 . Let 

2 r x 3 

a(x 2,6^3 ):=-- / c^i(x 2 ,C- 2 ,t)dt. 

1 Jo 

Since ci 1 is compactly supported, a is bounded, and one can form the unitary 
pseudo-differential operator exp(iU), A = Op™ (a). We have 

exp(-iA) Op™ (N h ) exp(z^) = Op™ (N h ) + exp(-iA)[Op™ (N h ) ,exp(L4)]. 

The symbol of [exp(— iA) Op™ (Nh) , exp(zAl)] is 

-e~ ia {N, e ia j + 0{h 2 ) = h{N h , a} + 0{h 2 ) = h{N 0 , a} + 0(h 2 ), 

1 

where Nq is the principal symbol of Nh, which satisfies: 

N 0 =e 3 +o(e 3 ). 

Therefore {Nn,a} = {d;o} + 0(d)- Since 

r o ^ ^ da 

{?3; a l = 2 C3 t^ = — ^3Cl,l, 

we get 

exp(-L4) Op™ (Nh) exp(L4) = Op%(JV h - h£ 3 c lyl + 0(h£) + 0(h 2 )), 

which shows that we can remove the coefficient of 3. The new operator 
given by the conjugation formula = exp(— iA) Op^ (Nh) exp(i^4) has a 
symbol of the form 

(4.8) JV' 11 = u 2 (x 2 , 6) (d + hxf) + hb(x 2 , C2, s(x 2 , 6)) + r h , 

where r h = 0(hx 3 ) + 0(h^ 3 ) + 0(d) + 0 (h 2 ). 

This proves Theorem 2.8. 

4.1.4. The second Birkhoff normal form. We now want to perform a Birkhoff 

r-i I u 

normal form for A [i relative to the “second harmonic oscillator” 

T 2 (x2-,&) (d + hx l) ■ 

Using Notation 2.10, we introduce the new semiclassical parameter h = h 2, 
and use the relation 

opsw: 1 *)=op^(n™-»). 

Thus, let d := The new symbol Nj^’^ has the form: 

(v_ 2 (x 2 , hi 2 )(H + ®§) + b(x 2 , /i£2, s(x 2 , h &)) 

+ h~ 2 rl 2 (x 2 ,hi 2 ,x 3 ,hi 3 ) S ) . 


N|d\^2,C2,®3,£3) = h 2 
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We introduce momentarily a new parameter fi and define 


Ni 1] ’ # (*2,|2,®3,l3;^) : = if 2 (^ 2 ,^ 2 )(II +x\) +b{x 2 ,ni 2 ,s(x 2 ,ni 2 )) 

+ h~ rl 2 (x 2 ,^2,x 3 ,h^ 3 ). 

Notice that £2, x 3 , £3; h) = h _2 N^’^(2;2, £2, x 3 , £3). We define now a 

space of functions suitable for the Birkhoff normal form in (x 3 , £3 ,h). Let us 
now use the notation of the Appendix introduced in (A.4) in the case when 
the family of smooth linear maps M 2 —y M 2 is given by 

^ t ,R2(^2,6) = (X2,ni2) ■ 


Let 

& := C(1) K 2, 

where the index M 2 means that we consider symbols on M 2 . More explicitly, 
we have 


& = {d s. t. 3c E 5(1; [0,1] x(0,1]) R 2 : d(x 2 , 6; H, h) = c(^ )R 2 (x 2 ,6); (i, h)} . 


Then we define 

£ ■= ^Ix3,i3,hj , 

endowed with the full Poisson bracket 


<§ x $ 3 (/, g) !->■ 


{f,d} = X] 

1=2,3 


df dg 
dij dxj 


dg_dj_ p 

dxj 


and the corresponding Moyal bracket [/, g ]. We remark that the formal Tay- 

MTU ~ 

lor series of the symbol N) ( ’ B (x2, &,x 3 , £3; //) with respect to (x 3 ,i 3 ,h) be¬ 
longs to S. We may apply the semiclassical Birkhoff normal form relative to 
the main term i/ 2 (x 2, + x\) exactly as in Section 3.2.1 (and also [25, 

Proposition 2.7]), where we use the fact that the function 

(x2,i2,x 3 ,i 3 -,n,h) (-A (^ 2 (x 2 ,//6)) _1 

belongs to $ because v} > C > 0 uniformly with respect to g. Let us 
consider 7 E $ the formal Taylor expansion of L -2 r^ 2 (x2, ^ 2 , x 3 , h£ 3 ) with 
respect to (x 3 , £3, h). The series 7 is of valuation 3 and we obtain two formal 
series k, t E § of valuation at least 3 such that 

[*. *3 + £3] = 0 


and 

e ^- lad ^ 2 (* 2 , ^2)(H + xl) + 7) = L 2 (.X2 , m6 )(H + X 2 ) + «. 

The coefficients of r are in 5(1) and one can find a smooth function 77 E 5(1) 
with compact support with respect to (x 3 ,£ 3 ,h) and whose Taylor series in 
(x 3 ,i 3 , h) is t. By the Borel summation, 77 will actually lie in 5(m / ) with 
m'(x 2 , £2) x 3 , £3) = {(x 3 , i 3 ))~ k for any k > 0, uniformly for small h > 0 and 
g E [0,1]. Notice that e C(m) with m = ((x 3 ,^ 3 )) 2 ^ 1, and that 

mm' = 0(1). 
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Then, we can apply Theorem A.3 with the family of endomorphisms of M 4 
defined 

7V,R 4 (^2,12,^3,13) = (x2,^i2,X3,i 3 ) ■ 

Thus, the new operator 

yjl h — e * h_1 °Pii Op” r h 

is a pseudo-differential operator whose Weyl symbol belongs to the class C (m) 
modulo h°°S( 1) (see the notations of Theorem 2.11). Moreover, thanks to 
Theorem A.4, its symbol M/j admits the following Taylor expansion (with 
respect to (£ 3 , £ 3 , h)) 

b( x 2 , s(x 2 , m6)) + T 2 (x 2 , m6)(H + x l ) + K - 

We write k = Ylm+ 2 t >3 c m,e{ x 2 , ^ 2 )\z 3 \* 2m h e and we may find a smooth 
function g*(x 2 , /i£ 2; ■Z', h) such that its Taylor series with respect to Z’, h is 

^2 c m ,e(x2,ix&)Z m h e . 

2ra+2£^3 

We may now replace g by h, which achieves the proof of Theorem 2.11. 

4.2. Spectral reduction to the second normal form. This section is 
devoted to the proof of Corollary 2.13. 

4.2.1. From to N In this section, we prove Corollary 2.9. 

Lemma 4 . 1 . We have 

n (a4 1]j , fan) = 0(h~ 2 ), n (a d iy ,= 0(h~ 2 ). 

Proof. The first estimate comes from Proposition 3.6 and Corollary 3.7. The 
second estimate can be obtained by the same method as in the proof of Corol¬ 
lary 3.7. □ 

Let us now summarize the microlocalization properties of the eigenfunc¬ 
tions of in the following proposition. 

Proposition 4 . 2 . Let xo be a smooth cutoff function on M being 0 in a 
neighborhood of 0 and let 6 E (0, ^). Let x be a smooth cutoff function 
being 0 on the bounded set {x| + &(x 2 , £ 2 , s(.x 2 , £ 2 )) A A)} and 1 on the set 

r-i 1 u 

{x 2 + &(x 2 ,£ 2 , s( x 2 , £ 2 )) A /?o + e}, with i > 0. If X is an eigenvalue of N_h 
such that A ^ foh and if if is an associated eigenfunction, then we have 

Ops (x (^, 6 , 3 : 3 ))^ = 

and 

o P £ (xo A = o(h°°)\\if\\. 

Proof. The proof follows exactly the same lines as for Propositions 3.11 and 
3.12. □ 
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Lemma 4.1 and Proposition 4.2 on the one hand and Propositions 3.11 and 
3.12 on the other hand are enough to deduce from Theorem 2.8 the point (a) 
in Corollary 2.9. The point (b) easily follows from Corollary 2.4. 

4.2.2. From to 971^. Let us now prove the point (a) in Corollary 2.13. 
We get the following rough estimate of the number of eigenvalues. 

Lemma 4.3. We have 

(4.9) N (mL 1],H , M 2 ) = N (St*, (3oh 2 ) = 0(h ~ A ), 

(4.10) n(q nl/3 0 h 2 ) =0(h~ 4 ). 

Proof. First, we notice that and are unitarily equivalent so that 

(4.9) holds. Then, given 77 > 0 and h small enough and up to shrinking 
the support of g * and by using the Calderon-Vaillancourt theorem (as in the 
proof of Lemma 3.4), ^ 971^ in the sense of quadratic forms, with 

= Op“ (h 2 b(x 2 , hi 2 , s(x 2 , hi 2 ))J + h 2 J h Op“ ((k 2 (x 2 , - gj . 

Since u 2 ^ c > 0, we get 

Op h ( h 2 b(x 2 ,hi 2 ,s(x 2 ,hi 2 ))) + /i 2 J4Op“ ((V(x 2 , ^ 6 )) - 17 ) 

^ Op (h 2 b(x 2 , hi 2 , s(x 2 , / 16 ))) + |^ 2 Jh • 

We deduce the upper bound (4.10) by separation of variables and the min- 
max principle. □ 

The following proposition deals with the microlocal properties of the eigen¬ 
functions of 

Proposition 4.4. Let rj G (0,1), 5 G (0, ,C > 0. Let x be a smooth 
cutoff function being 0 on {b(x 2 ,£ t2 ,s(x 2 ,f 2 )) ^ /?o} and being 1 on the set 
{b(x 2 ,£ t2 ,s(x 2 ,£ i2 )) ^ /3q + i}, with e > 0. Let also xi be a smooth cutoff 
function on M 2 , being 0 in a neighborhood of 0 . 

If A is an eigenvalue offfsuch that A ^ /3q h? and iff) is an associated 
eigenfunction, we have 

(4.11) Op£ (x(x2, /*&)) if = 0(h°°)\\if\\ 

r-i I u 

and if X is an eigenvalue o/92j t such that A ^ boh 2 + Ch 2+ri and if if is an 
associated eigenfunction, we have 

Op'h ( Xi(h~ 5 (x 3 ,i 3 ))) 


(4.12) 


if = 0(h°°)\\if\\. 
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Proof. The estimate (4.11) is a consequence of Proposition 4.2. Then, let us 

r-i I u 

write the symbol of Tt) ( : 

#' = h 2 v 2 (x 2 ,h£ 2 ) + x 2 ^j+h 2 b{x 2 ,hf 2 , s{x 2 ,hf 2 ))+li h2 (x 2 ,hf 2 ,x 3 ,hf 3 ). 

We write 


^h U °Ph(xi(h S (x 3 ,£ 3 ))\ if,Op% (xi{h 5 (x 3 


= A ll OPh (Xi{h \x 3 , 


+ 




'\OVh[xi{h 5 (x 3 ,£ 3 ))) ,OVh(xi{h S {x 3 , 


We get 



y±h U ,OPh (xi(h 5 {x 3 , 


,Op^(xi(/t \x 3 , 


if 


^ Ch 3 


Op h (X^h 5 (x 3 ,&))) if 


where we have used (4.11). Then, we use that 

b(x 2 , h£ 2 , s(x 2 , h £ 2 )) ^ bo, n 2 (x 2 , h£ 2 ) ^ c 0 > 0, A ^ b 0 h 2 + Ch 2+V , 
and the Garding inequality to deduce 


h 2 ( Ch 25 - Ch v 


Op h [xi(h 5 (x 3 ,£ 3 ))) if 
^ Ch, 3 


Op)C (x 3 , 

The desired estimate follows by an iteration argument. 




□ 


In the same way we can deal with 

Proposition 4.5. Let r] E (0,1), 5 E (0, , C > 0. Let x be a smooth 

cutoff function being 0 on {b{x 2 ,f 2 ,s(x 2 ,f 2 )) ^ /3o} and, being 1 on the set 
{b(x 2 , £ 2 , s(x 2 , £ 2 )) ^ fto + e}, withe >0. //A is an eigenvalue ofD)such 
that A ^ /?o^ 2 an d if if is an associated eigenfunction, we have 

(4.13) Op^ (x(^2, ^£ 2 )) ^ = 0(h°°) Ill'll 

and if A is an eigenvalue of 93 tjj such that A ^ boh 2 + Ch 2+ri and if if is an 
associated eigenfunction, we have 

(4.14) Op£ (xi(h- 5 (x 3 ,£ 3 ))) if = 0(h°°)\\if\\. 

Proof. In order to get (4.13), it is enough to go back to the representation 
with semiclassical h, that is 971J = A4|. Indeed the microlocal estimate 
follows by the same arguments as in Propositions 3.11 and 3.12. Then, (4.14) 
follows as in Proposition 4.4. □ 
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Propositions 4.4 and 4.5 and Theorem 2.11 standardly imply the point (a) 
in Corollary 2.13. 

4.2.3. From to Let us now prove the point (b) in Corollary 2.13. 

Note that the point (c) is just a reformulation of (b). 

Let us consider the Hilbertian decomposition where 

the symbol IVlJ^ of is 

h 2 b(x 2 , h£ 2 , s(x 2 , h£ 2 )) + (2k - 1 )h 3 v 2 (x 2 , hf 2 ) + h 2 g*(h, (2k - 1 )h,x 2 , h| 2 ). 
There exists ho > 0 such that for all k ^ 1 and h £ (0, ho), 

(M [ * U iP,iP) 

> <°Pfe (ji 2 b(x 2 ,h£ 2 ,s(x 2 ,h£ 2 )) + (2k - l)h 3 (v 2 (x 2 ,h£ 2 ) - s)') ip,ip). 

\k] W 

Since each eigenfunction of associated with an eigenvalue less than 

/3oh 2 provides an eigenfunction of 9JT^, we infer that the eigenfunctions of 
are uniformly microlocalized in a (x 2 , £ 2 )-neighborhood of (0,0) as 
small as we want. Therefore, on the range of 1(_ 00 ^/^(DJtJ^), we have 

(m [ £ u ip,ip) 

> (Op™ (h 2 6(x 2 ,h£ 2 ,s(x 2 ,h£ 2 )) + (2k - l)h 3 (z/ 2 (0, 0) - 2e)) ip,ip). 
and, with the Garding inequality in the h-quantization, we get 

(9J lf^ip,ip) ^ (Op^ (h 2 fe 0 + (2k - l)h 3 (z/ 2 (0,0) - e) - Ch A ) ip, ip). 

This implies the point (b) in Corollary 2.13. 

5. Third Birkhoff normal form 

5.1. Birkhoff analysis of the first level. In this section we prove Theo¬ 
rem 2.15. 

We consider = Op™ with 

M n ],i = hb(x 2 ,&,s(x 2 ,&)) + hi^ 2 (x 2 ,6) + hg*(tfi,h5,x 2 ,&)- 

By using a Taylor expansion, we get, 

(5.1) 

M^ A = hb 0 + ^Hess( 0)0 )6(x 2 ,^ 2 ,s(x 2 ,^ 2 )) + h5^ 2 (0,0) + cx 2 hi + d£ 2 hs 

+ hO((hKz 2 ) 3 ), 

where c = d X2 v 2 (0, 0) and d = d^ 2 v 2 ( 0,0), and we have identihed the Hessian 
with its quadratic form in (x 2 ,£ 2 ). 
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Then, there exists a linear symplectic change of variables that diagonalizes 
the Hessian, so that, if L% is the associated unitary transform, 


L1M 


L'r# = Op!" 


M, 


[i]J 


with 

= hbo + 7 ^ 0 (x 2 + il) + 0,0) + cx 2 h 5 + d£ 2 hi + hO((h^, z 2 ) 3 ), 

where 

0 = yj det Hess (0)0 )6(x 2 ,6, s(x 2 ,&)) . 

Since (d X3 b)(x 2 , g 2 , s(x 2 , £ 2 )) = 0 and (0,0) is a critical point of s, we notice 
that d% 2X3 b( 0,0,0) = df 2a . 3 &(0,0,0) = 0. Thus 

det Hess (0i0 ,o)&(0, 0, 0) = 0 2 d 2 3 b(0,0 ,0). 

Using that b is identified with box (see Remarks 2.2 and 3.1), this provides 
the expression given in (1.14). 

Note that c 2 + d 2 = || (V X2 ^ 2 ^ 2 )(0, 0)11 2 since the symplectic transform is 
in fact a rotation. Moreover we have 


0{x 2 ~\-^, 2 )~^rcx 2 fi 2 +d^- 2 h 2 




dh ] A 2 \ & + & 

-t)) - fi — 


Thus, there exists a unitary transform U i, which is in fact an h-Fourier 

Integral Operator whose phase admits a Taylor expansion in powers of h 2 . 
such that 

U\LlM™L h U h i =: F h = Op£ (F h ), 

where 


F h = hb 0 + H 2 v 2 (0, 0) - 


26 


h 2 + h( w -\z 2 \ 2 + 0 {{h~ 2 ,z 2 f' 


Now we perform a semiclassical Birkhoff normal form in the space of formal 
series M[ai 2 , £ 2 , equipped with the degree such that x\fJi'h 2 is t + m + n 
and endowed with the Moyal product. Let Fj be the full Taylor series of F_n- 
We hnd a formal series t(x 2 , ^ 2 > ) with a valuation at least 3 such that 


Ah x ad 7 


Fl = F£, 


where F^ is a formal series of the form 


Fj = Hbo + ft 2 iZ (0, 0) 



22 


and k T is a formal series in M[h 2 , \z 2 \ 2 } (and that can be also written as a 
formal series in Moyal power of \z 2 \ 2 , say ( k T )*). 
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Let f(x2, ^2) m) be a compactly supported function whose Taylor expansion 
at (0,0,0) is equal to t{x 2, £ 2 , A 4 )- By the Egorov theorem A.2, uniformly with 
respect to the parameter p , we obtain that 

e -m-iOpJ»(f) Op^(F /i 2 )e a " 1 ° p ^ (f) =: Op^(F M ) 

is an ^,-pseudo-differential operator depending smoothly on p. Expanding 
Ffj in powers of p in the 5(1) topology, and letting /i = y/h, we see that 
^Vh = F h + Gh, where 


F h = hb 0 + /i.21/ 2 (0, 0) 



with k a smooth function with a support as small as desired w.r.t. its second 
variable, and Gh = hO(\z 2 \°°). It remains to notice that Op™ [k{hz , l^) 2 )^ 

can be written as k*(h 2 , Kf) modulo Op™ (0(\z2\°°)). This achieves the proof 
of Theorem 2.15. 


5.2. Spectral reduction to the third normal form. Corollary 2.16 is a 
consequence of the following lemma and proposition. 

Lemma 5.1. We have 

n ( M™,m ) = o(tr 2 ), n (F n ,b 0 h + Ch 1+ v) = o(/r 1+r? ). 

Proof. The first estimate follows from Lemma 4.3 and the second one from a 
comparison with the harmonic oscillator in x' 2 - □ 

The last proposition concerns the microlocalization of the eigenfunctions. 

Proposition 5.2. Let 77 E (0,1), 6 E (0, , C > 0. Let x be a smooth cutoff 

function being 0 in a bounded neighborhood of (0, 0) and 1 outside a bounded 
neighborhood o/(0,0). If A is an eigenvalue of or of Fn such that 

A ^ &o^ + Ch 1+V and if if is an associated eigenfunction, we have 

OpJ? (x(/r*(x 2 ^2)))it> = o{h°°). 

Proof. The proof is similar with the one of Proposition 4.4. □ 


Appendix A. Egorov theorems 

We start with the classical result (see for instance [31, Theorem 11.1] 
and [26, Theoreme IV.10]). 

Theorem A.l ([31, Theorem 11.1, Remark (ii)]). Let P and Q be h-pseudo- 
differential operators on M d , with P E Op™ (5(1)) and Q E Op™ (5(1)). 
Then the operator e?A Pe~h® is a pseudo-differential operator in Op)]' (5(1)), 
and 

esAPe - ^ — Op™(p° k) E /iOp™ (5(1)) . 

Here p is the Weyl symbol of P, and the canonical transformation n is the 
time-1 Hamiltonian flow associated with principal symbol of Q. 
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From this classical version of Egorov’s theorem, one can deduce the fol¬ 
lowing refinement that is useful when p does not belong to £(1) (as it is the 
case in this paper). 

Theorem A. 2. Let P and Q be h-pseudo-differential operators on with 
P E Op)) (S(m)) and Q E Op))) ( S(m')), where m and m! are order functions 
such that: 

(A.l) m 1 = 0(1); mm! = 0(1). 

Then the operator Pe~hQ is a pseudo-differential whose symbol is in 
S(m), and e^Pe~^ — Op f(p o k) E /i Op)) (5( 1)). 

Proof. The proof is based on the following observation. In order to compare 
Op f(p o k!) and e^Pe - ^, we consider the derivative: 

^ [e^ Q Op))(po K t ~ T )e~ ! ^ Q ^ 

= e £0 ^ [Q, Op)) (p o n t ~ T )\ + ± Op))(po K ^)\ e~^. 

From Hypothesis (A.l), the term [Q, Op)) (p o n t ~ T )] belongs to Op)) (5(1)); 
moreover, if we denote by qo the principal symbol of Q, we have 

^ Op))(po k*-t) = - Oph({qo,P o 

which implies that this term is also in Op)) (5(1)). By symbolic calculus, we 
see that 

(A.2) l[Q , Op f(p o rJ~ T )\ + ‘{ r Op,J)(p o «‘- T ) E h Op)) (5(1)), 

uniformly for t. t in compact sets. It follows by integration from 0 to t that 

(A.3) e^Pe - ^® = Op))(p o if) + h f Pi(s)e~^^ds, 

Jo 

for some Pi(s) E Op)) (5(1)), uniformly for s E [0, t]. Applying Theorem A.l 
to the integrand, we see that e^Pe”^ — Op)) (p o n!) E hOp)) (5(1)). □ 

In order to quantize the formal Birkhoff procedure of Section 4.1.4, one 
needs to consider symbols in a class C stable under the Moyal product. For 
that purpose we first define the families of symbols 5(m; [0,1] x (0,1]), that 
is of smooth functions a : M 2d x [0,1] x (0,1] —> C such that, for any a E N 2d , 
there exists C a such that, V(z; p, h) E R 2d X [0,1] x (0,1], 

\dfa(z] p, h)\ ^ C a m(z) 

and where m is an order function on M 2<i . The pair (p, h) is considered as a 
parameter. 
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Then, let (<^/i)^e[o,i] be a smooth family of linear maps M 2d —> M 2d and 
define the following families of symbols on M 2rf by 

(A.4) C(m) = |a G 5(m; [0,1] x (0,1]); a(z; p,h) = afo^z); p,h) 

with a G S(m; [0,1] x (0,1])| . 

Theorem A.3. Let P and Q be h-pseudo-differential operators on with 
P G Op^ ( C(m )) and Q G Op)]] ( C(m')), where m and m' are order functions 
such that: 

m ^ 1; m' = 0(1); mm' = 0(1). 

Then e^Pe~^ = P + R, where P G Op)) ( C(m)), R G h°° Op)) (5(1)), 
and with P — Op))(po n) G hOp)) (C(l)). 

Proof. Since is linear, one can see (using for instance [31, Theorem 4.17]) 
that C is stable under the formal Moyal product, i.e. for all order functions 
m\ and m 2 , we have 

(i C(m\)) * ( C(m 2 )) C C(m, 1 m 2 ) + h°°S(l). 

Let k be the canonical transformation associated with Q. Then, since 
m ^ 1, we have jjokG C(m); indeed, if we write the Hamiltonian flow of Q 
in terms of the variable 5 = (pffz), we see from the linearity of that the 
components of the transformed vector field belong to C(m'). Therefore ok 
is of the form k^ o for some diffeomorphism k^ depending smoothly on p. 

Therefore, both terms in (A.2) belong to Op)) (C(l)). Applying this argu¬ 
ment inductively in (A.3), we may write, for any k > 0, 

e iQp e -iQ - Op J)(p ok)- (hPi + h 2 P 2 + • • • + h k P k ) G h k+1 Op)) (5(1)), 

with Pj G Op)) (C(l)). By a Borel summation in h, parametrized by z = 
Pfi(z), we can find a symbol P G Op)) (C(l)) such that we have the asymp¬ 
totic expansion in Op)]' (5(1)): 

P ~ hP\ + h 2 P 2 + • • • 

We conclude by letting P = Op)]’ (po k) + P. □ 

We will also need to examine how the Egorov theorem behaves with re¬ 
spect to taking formal power series of symbols. For this, it is convenient to 
introduce a filtration of S(m). 

Theorem A.4. Let m be an order function on M 2d ; and let (Oj)j^ be a 
filtration of S(m), i.e.: 

Oo = S(m), O j+ 1 C Oj. 

Let P = Op)] 1 p and Q = Op)]’ q be h-pseudo-differential operators on M d , 
with p G S(m) and q G S(m'), where m' is an order function such that m' 
and mm' are bounded. 
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Assume that: 

(A.5) ^ &d q (Oj) C Oj + i] Vj^O. 

T/ien for any k ^ 0, i/ie Wet// symbol of the pseudo-differential operator 
e^Pe~^ — Yl k j=o j\(h a dQ) J -P belongs to Op^ (CV|_i). /n oi/ier words, the 
series of ex ad q)P converges to Pe~hQ for the filtration (Oj)j & ^. 

Proof. By the Taylor formula, we can write 

e^ Q Pe-t Q = ^ ^(&d ih -iQ) j P+^{ad ih -i Q ) k+1 f (1 -t) k e^ Q Pe~^ Q dt. 
j= o j! Jo 

By Theorem A.2, we see that the integral belongs to Op^ ( S(m )) = Op“ (Oq). 
Therefore, by Assumption (A.5), the remainder in the Taylor formula lies in 
Op w h (O k+1 ). □ 
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